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Abstract 

Assuming gauge-mediated supersymmetry (SUSY) breaking (GMSB), we simulate precision 
measurements of fundamental parameters at a 500 GeV e + e~ linear collider (LC) in the scenario 
where a neutralino is the next-to-lightest supersymmetric particle (NLSP). Information on the 
SUSY breaking and the messenger sectors of the theory is extracted from realistic fits to the 
measured mass spectrum of the Minimal SUSY Model (MSSM) particles and the NLSP lifetime. 



Contribution to the Workshops: 
2 nd ECFA/DESY Study on Physics and Detectors for a Linear e + e~ Collider 
Lund, Frascati, Oxford, Obernai - June 1998 to October 1999 



Supersymmetry must be broken if it is to describe nature, and GMSB [1] is one attrac- 
tive way to realize this, also providing natural suppression of the SUSY contributions to 
flavour- changing neutral currents at low energies. In GMSB models, the gravitino G is the 

LSP with mass given by txiq — — 2.37 ( iocrr e v ) wnere * s the fundamental 

SUSY breaking scale. The GMSB phenomenology is characterised by decays of the NLSP 
to its Standard Model partner and the G with a non-negligible or even macroscopic life- 
time. In the simplest GMSB realizations, depending on the parameters M mess , N mess , A, 
tan/?, sign(/i) defining the model, the NLSP can be either the lightest neutralino Ni or 
the light stau f\. For this study [2], we generated several thousand GMSB models fol- 
lowing the standard phenomenological approach [3] and focused on the neutralino NLSP 
scenario, for which we selected several representative points for simulation. Our aim was 
to explore the potential of a LC in extracting the fundamental model parameters. 

Firstly, we investigated the sensitivity in determining the GMSB parameters at the 
messenger and electroweak scales from the knowledge of the sparticle masses that could 
be obtained from threshold-scanning techniques. We used a sample model with 
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Parameter 


Fitted value 


-^mcss 


(161 ± 2) TeV 


A 


(76.01 ± 0.08) TeV 


1 ' mess 


0.9994 ±0.0009 


tan (3 


3.50 ±0.03 



Table 1: Results of fits to the parameters of the GMSB model described in the text, 
starting from a possible set of light sparticle masses measurements from threshold scans. 
A 200 fb' 1 run at the LC is assumed. 

M mess = 161 TeV; N mess = 1; A = 76 TeV; tan/? = 3.5; fi > 0, producing a rather light 
sparticle spectrum, and assumed a total of 200 fb" 1 collected between 200 and 500 GeV 
c.o.m. energies at a LC. By just considering the shape of the total cross sections for 
several kinematically allowed SUSY production processes as functions of y/s close to the 
thresholds, we inferred the following approximate precisions for the sparticle masses: 



A(m^) ~ 0.2 GeV; A(m^ 2 ) ~ 0.8 GeV; A(m dl ) ~ 0.1 GeV; 

A(m e ~J ~ 0.2 GeV; A(m e ~J ~ 0.2 GeV; A(m M J ~ 0.8 GeV; (1) 

A(m fl ) ~ 0.8 GeV; A(m f2 ) ~ 2.0 GeV; A(m h o) ~ 0.1 GeV. 

By performing a fit to minimise a x 2 based on these errors, with the true (model- 
dependent) values as the central ones in the fit, we obtained an estimate of the precisions 
on the underlying parameters, as shown in Tab. 1. We checked that these precisions are 
typical for the class of models we considered. 

Then, we considered N\ lifetime measurements in the whole allowed cr^ range, per- 
forming event simulation in detail for our set of representative GMSB models. Indeed, 
since the N\ lifetime is related to \fF by 

_16nVF\_ l( VF \ 4 / 
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the GMSB framework provides an opportunity to extract information on the SUSY 
breaking sector of the theory from collider experiments that is not available, e.g., in 
supergravity-inspired models. 

Typical neutralino lifetimes for our models range from microns to tens of metres. While 
the lower bound on cr^ i comes from requiring perturbativity up to the grand unification 
scale [3] , the upper bound is only valid if the G mass is restricted to be lighter than about 
1 keV, as suggested by some cosmological arguments [4] (cfr. Fig. 1). 

For given Ni mass and lifetime, the residual theoretical uncertainty on determining 
\[F is due to the factor of order unity B in Eq. (2), whose variation is quite limited in 
GMSB models (cfr. Fig. 2a). 

In addition to the dominant Ni — > 7G decay, it was fundamental to our analysis 
to take the A r 1 — > Gff decays into account, in order to use the tracking detectors for 
measurements of shorter neutralino lifetimes. We performed a complete study of these 
channels and found that in most cases of interest for our study the total width is given 
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Figure 1: Scatter plot of the neutralino NLSP lifetime as a function of the messenger scale 
M mess (a) and (b). For each set of GMSB model input parameters, we plot the lower 
limit on ct^ , corresponding to \[F ~ v^Fmess — \MM mcss . We use only models that fulfil 
the limit uiq < 1 keV =>- VF mcss ^ y/~F < 2000 TeV suggested by simple cosmology. 

approximately by 

+ r(N 1 ^ZG)B(Z^ff), 

(3) 

where the expressions for the widths of the 2-body N\ decays are well-known [5]. In 
Fig. 2b, the branching ratio (BR) of the Ni — > 7G decay is compared to those of iVi — > ZG 
and N\ — > h°G (in the on-shell approximation) and those of the main Ni — > ff'G channels 
(including virtual-photon exchange contributions only). 

To generate GMSB events, we modified SUSYGEN 2.2/03 [6], to take the 3-body neu- 
tralino decays into account as follows. We implemented in CompHEP 3.3.18 [7] a home- 
made lagrangian including the relevant gravitino interaction vertices in a suitable approx- 
imation. Then, we studied the kinematical distributions of the Ni — > ffG channels and 
passed the results to the event generator numerically. For each sample GMSB model, we 
considered in most cases a LC run at a c.o.m. energy such that the only SUSY production 
process open is NLSP pair production e + e~ — > N1N1, followed by iV\ decays through all 
possible channels. For more challenging models where the light SUSY thresholds are close 
to each other, we simulated also events from i?-slepton pair production and used some 
selection to isolate the N1N1 events, for which the Ni production energy is fixed by the 
beam energy (we also took into account initial-state radiation as well as beamstrahlung 
effects), allowing a cleaner cr^ measurement. 

The primary vertex of the events was first smeared according to the assumed beamspot 
size of 5 nm in y, 500 nm in x and 400 /iin in z and then the events were passed through a 
full GEANT 3.21 [8] simulation of the detector as described in the ECFA/DESY CDR [9]. 
The tracking detector components essential to our analysis included a 5-layer vertex de- 
tector with a point precision of 3.5 /zm in r<f> and z, a TPC possessing 118 padrows with 
point resolution of 160 fim in r<fr and 0.1 cm in z. In addition, we assumed an electromag- 
netic calorimeter with energy resolution given by (10.3/\^E + 0.6)%, angular pointing 



T(N, - ffG) ~ T(iV 1 - 7 G)^A7Q 
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Figure 2: (a) Scatter plot showing the relation between the neutralino NLSP lifetime and 
the fundamental scale of SUSY breaking v~F, i.e. the factor B in Eq.(2), as a function of 
the neutralino mass in GMSB models of interest for the LC (100 < < 250 GeV). Big 
grey dots in represent models with a light R-selectron (102-150 GeV), small black dots 
are for the heavier selectron case (150-430 GeV). (b) Scatter plot for the BR's of various 
N\ decay channels as a function of the N\ mass. Dots in different grey scale (colours ) refer 
to the decays Ni — > jG, jVj — > ZG (including off-shell effects), and to hadrons or e + e~ 
plus gravitino via virtual photon, as labelled. For reference, we also report results for the 
2-body Ni — > h°G decay in the on-shell approximation, whose BR is always negligible. 

resolution of 50/\/E~ mrad and timing resolution of 2/\^E ns. The dimensions of the 
whole calorimeter (electromagnetic and hadronic) were 172 cm < r < 210 cm and 280 cm 
< \z\ < 330 cm. 

The probability for a single neutralino produced with energy Eft to decay before 
travelling a distance A is given by P(X) = 1 — exp(— X/L), where L = CTft l (f3j)ft l is the 
Ni "average" decay length and ([3 r y)ft l = (Eft. /m^ — l) 1 / 2 . 

For L less than a few cm, we used tracking for measuring the vertex of N\ — > Gff 
decays. When L is very short, less than a few hundred /zm, the beamspot size becomes 
important and a 3D procedure is not appropriate. Instead, the reconstructed vertex was 
projected onto the xy plane, where the beamspot size is very small, and we used the 
resulting distributions to measure the Ni lifetime. We studied several GMSB models 
with crft in the allowed range and found that the intrinsic resolution of the method was 
approximately 10 /im. An example of the reconstructed 2D decay length distribution for 
a challenging model where the neutralino lifetime can be very short [2] is shown in Fig. 3 
for statistics corresponding to 200 fb _1 (r is the xy component of A). 

For 500 yum < L < 15 cm, we used 3D vertexing to determine the decay length 
distribution and hence the lifetime of the N\. Vertices arising from Ni — ► 7G and photon 
conversions in detector material were essentially eliminated using cuts on the invariant 
mass of the daughter pairs together with geometrical projection cuts involving the mass 
of the Ni and the topology of the daughter tracks. Methods of measuring the Ni mass 
using the endpoints of photon energies or threshold techniques, together with details of 
the projection cuts have been described [2]. Using 200 fb _1 of data, we concluded that a 




Figure 3: (a) Reconstructed projective radial distances, r, of the Ni — > ffG decay vertex, 
for a LC run on a short-lifetime model with L = 10 fim together with a fit to an exponential 
plus constant, (b) Shows the results of the corresponding fits for a range of L values. 



CTft measurement with statistical error of ~ 4% could be made using this method. 

For L larger than a few cm, we used the N% — > 7G channel, providing much larger 
statistics. The calorimeter was assumed to have pointing capability, using the shower 
shapes together with appropriate use of pre-shower detectors. Assuming the pointing an- 
gular resolution mentioned above, we demonstrated [2] how a decay length measurement 
can be made. We concluded that for lifetimes ranging from approximately 5 cm to ap- 
proximately 2 m this method worked excellently, with statistical precisions ranging from a 
few % at the shorter end to about 6% at the upper end of the range. We also investigated 
the use of timing information to provide a lifetime measurement, but found it to be less 
useful than calorimeter pointing. However, the use of timing might be relevant to assign 
purely photonic events to bunch crossings and to reject cosmic backgrounds. 

For very long lifetimes, we employed a statistical technique where the ratio of the 
number of one photon events in the detector to the number of two photon events was 
determined as a function of CTfj . This allowed a largely model- independent measurement 
out to crfj ~ few 10's m. The possibility of using the ratio of the number of no-photon 
events to one photon events was also discussed [2]. The latter allows a greater length 
reach, but relies on model-dependent assumptions. 

In Fig. 4, we summarise the techniques we have used as a function of L for a sample 
model. The criterion for indicating a method as successful is a measurement of L and 
the Ni lifetime to 10% or better. It can be seen that L can be well measured for 10 's of 
/im < L < 10's of m, which is in most cases enough to cover the wide range allowed by 
theory and suggested by cosmology. 

With reference to Eq. (2), we note that a 10% error in cTft corresponds to a 3% error 
in yF. This is of the same order of magnitude as the uncertainty on the factor £>, which 
parameterises mainly the different possible Ni physical compositions in GMSB models 
(cfr. Fig. 2a). We also checked explicitly that, in comparison, the contributing error 
from a neutralino mass measurement using threshold-scanning techniques or end-point 
methods is negligible [2]. 
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Figure 4: Summary of the techniques used here for a cr^ measurement at the level of 
10% or better. 

Hence we conclude that, for the models considered and under conservative assump- 
tions, a determination of \fF with a precision of approximately 5% is achievable at a LC 
by only performing N\ lifetime and mass measurements in the context of GMSB with 
neutralino NLSP. Less model dependent and more precise results can be obtained by 
adding information on the N\ physical composition from other observables, such as N\ 
decay BR's, cross sections and distributions. 
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